Current-spin coupling for ferromagnetic domain walls in fine wires 
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The coupling between a current and a domain wall is examined. In the presence of a finite current 
and the absence of a potential which breaks the translational symmetry, there is a perfect transfer 
of angular momentum from the conduction electrons to the wall. As a result, the ground state is 
in uniform motion. This remains the case when relaxation is accounted for. This is described by, 
appropriately modified, Landau-Lifshitz-Gilbert equations. 

PACS numbers: 75.60.-d, 73.40.Cg, 73.50.Bk, 75.70.i 



Spintronic devices have great technological promise 
but represent a challenging problem at both an applied 
and fundamental level. It has been shown theoretically|l|, 
2] that the direction of a magnetic domain might be 
switched using currents alone. Devices designed to use 
this principle often consist of multi-layers of magnetic 
and non-magnetic conductors. The advantages of similar 
devices based upon the current induced displacement of 
a domain wall are simplicity and the fact that the switch- 
ing current is much smaller Experimentally the 
current induced displacement of a domain wall has been 
clearly demonstrated and in recent experiments^,!^ the 
velocity of the wall was measured. 

The current induced motion of a magnetic domain in- 
volves the transfer of angular momentum from the con- 
duction electrons. The early theory 0,0 and most of the 
subsequent work0 is based upon some type of assump- 
tion about this torque transfer process and there is no 
real consensus on how this should appear in the (Landau- 
Lifshitz-Gilbert) equations of motion 0- The purpose of 
this Letter is to develop a complete first principles the- 
ory of this process for a domain wall, this based upon a 
specific model Hamiltonian and physically justified ap- 
proximations. This is intended to form the basis for the 
more complicated problem of multi-layers. 

Although similar conclusions are valid for the Stoner, 
and related, models, here attention will be focused upon 
the s — d-exchange model. The direction of the local mo- 
ments, Si, which make up the domain wall are specified 
by the usual Euler angles 9i and 4>i ■ To make diagonal the 
interaction +JSi ■ s*j, at site i, the axis of quantization of 
the conduction electrons, Sj, is rotated along this same 
direction. If ij}(fi,t) is the conduction electron spinor 
field then this amounts to making a SU(2) gauge trans- 
formation, if)(n,t) -> r(6i,<f>i)ip(fi,t), where r(9 i ,4> l ) = 
e #, Sje ,8, %e ,^ Sz = ( cos (0./2) + isin^/^sin^o-s - 
i sin(0j/2) cos 4>iO~ y ) and where a are the Pauli matrices. 
If there was SU(2) gauge invariance, there would exist 
three gauge particles, the analogy of the W ± and Z which 
mediate the weak interaction. However here, without 
this invariance, this transformation still introduces three 
gauge fields. The longitudinal such field has been ex- 
ploited in the development of theories of the Hall effect ||. 



It is Bazaliy et al. , who first developed a theory of the 
torque transfer to a magnetic domain, for the half-metal 
limit, based upon this same field. Their formalism gener- 
ates a transverse field in the Landau-Litshiftz equations 
through a finite derivative dV{9i, <f>i}/d<f>i where V{9i, 
is the expectation value of the energy as a function of the 
angles {9i,<j>i}. However, a finite ground state expecta- 
tion value for such a derivative implies that the solution 
is not stable and it follows that their finite velocity solu- 
tion must relax to one which is stationary with a finite 
common value of the 4n = (f>o, this anticipating the con- 
clusions of Tatar a and Kohno[Io|. 

Here it will be shown that the torque transfer effects 
can be correctly accounted for in a simpler U(l) theory in 
which the only generators are Si y for rotations about the 
axis perpendicular to the wall plane. The local spins are 
usually co-planer, see below, and the rotations r(9i) = 
e i9is v _ ( cos (#./2) — ism(9i/2)o- y ) are all that is needed 
to diagonalize +JS ■ s. This simpler approach can only 
generate a transverse gauge field, i.e., the one ignored in 
the earlier workj^, and which alone is found to be the 
origin of the torque transfer process. The longitudinal 
field, along with dV{9i, <j>i}/d<j>i, and the ground state 4>o 
are strictly null, and the correctly relaxed ground state 
does have a finite velocity in the absence of pinning. 

The s — d-exchange Hamiltonian is, 
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where c itJ is the conduction electron creation operator 
for spin a and site i. The uniform hopping integral is 
Ujaa' — tS aa > (the strange notation is for later use), and 
N is number operator. The exchange constants are J and 
J° for the coupling between the local moments and the 
conduction electrons and for the direct coupling between 
local moments. With the anisotropy constant A > 
the z direction is an easy axis while, with K\ > 0, the 
y axis is hard. These anisotropy constants include both 
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intrinsic anisotropy plus the effects of a demagnetizing 
field. 

It is first necessary to describe the domain wall. The 
Holstein-Primakoff transformation^]] 5 iz = 5 — b\b^ 

St = (25 - &i& i ) 1/2 6 i « (2S') 1 / 2 6 l is used to quan- 
tize the spin degrees of freedom using the direction 
defined by 9i and <f>i as the axis of quantization for 
the local spins. To within a constant, for the spin 
Hamiltonian in Eqn. Q, V{9i,<pi} = -AS 2 cos 2 9 t + 
K^S 2 sin 2 9i sin 2 & — 2<y> cos where ^4 = [(25 — 
1)/25]A°, if_L = [(25 - l)/25]/^, and % is the an- 
gle between the spins at sites i and j. The energy 
K±S 2 sin 2 9 i sin 2 (f>i favors a wall with fa = 0, whence 
assuming a slow variation in the angle 9i, i.e., that 
9ij w V^i • r*y , where = r*j — is the vector which 
joins near neighbor planes of the wall in the z direc- 
tion (fij = a), and again dropping a constant V{#i} = 
+ J s 5 2 (V6»j • fij) 2 - AS 2 cos 2 9i. The domain wall struc- 
ture, which lies in the x — z-plane, can be determined 
by minimizing this V{#;}, and the standard result is 
9{zi) — 2cot _1 e~( Zi ' w \ where Zj is coordinate of site 
i. The wall width w = a (Js/2^1) 1 / 2 . The local spins 
rotate so 9 = becomes 9 = tt with increasing z. 

In order to have precise expressions for the conju- 
gate position and momentum of the wall, it is ob- 
served that the later is the generator of displace- 
ments. The product of small y-axis rotations R(Az) = 
Iliexpm - 0j)(Az/a)(S iy /h)) w (1 + (1/2)£.(V0 4 • 
fij)A(z/a)(2S) 1 ^ 2 (bl — bA), produces a wall translation 
by Az, and so the momentum, per atom in a plane of the 
wall, p z = i(h/2)(v c /a 2 A)j: i (V9 l ■ r u )(25) 1 / 2 (6j - bA, 
where v c is the cell volume and A is the cross-sectional 
area of the wire. The Goldstone boson which restores the 
translational invariancc of TL can then be identified as, 

b\ = {J/BA)V\v c /aA) ^(Vfi, • fj) b\, 

i 

where the constant of proportionality reflects the require- 
ment that [b w ,b\^\ — 1 and uses the explicit wall so- 
lution for 6{z). (The modes which disperse from this 
Goldstone boson are excitations of the wall.) Thus 
p z = ih{2A/ J) 1 l i {2S) 1 / 2 {l/a){bl - b w ) while the con- 
jugate coordinate, such that [i^Pz] = Hi, is zq = 
a(J/2A)Vy/2)-V 2 (bi + bJ. 

Returning to the description of the conduction elec- 
trons, performing the above described SU(2) rotations, 
the conduction electron part of Eqn. reduces to, 

Ti e = - 51 (clo-tijvv'Cjv' + H - c -) ~ JS X! Siz ~ ^N, 

<ij>aa' i 

(2) 

in which now tij — tr l (9i,4>i)r{9j,(j)j). But since 4n = 
this simplifies to the U(l) result: 

Uj = t(cos(9 l /2)+isixi{9 i /2)(T y )(cos(9 j /2)-isixi(9.j/2)a y ). 

(3) 



which is, Uj = t[cos(6»i/2) cos(^-/2)+sin(^/2) sin(0,-/2)+ 
i(sin(9i/2) cos(9j/2) — sin(8j/2) cos(9i/2)a y ], using cr y 2 = 
1. The standard identities then lead to fy = t[cosV9i ■ 
+ i(S79i ■ fij)a y ], and correct to second order in the 
gradient, this can be written as, 

Uj = e ^ A ' d? t cos( V0i -fij); A= (V0< • fij)s y z (4) 

In order that this central result be valid, the conduc- 
tion electrons must follow the local spin magnetization 
as they pass through the wall and this requires that the 
adiabatic theorem be satisfied. When the transverse field 
represented by Atj^ — it(V9i-fij)cr y is ignored (the equiv- 
alent to A = 0), the eigenstates of TL e are also eigenstates 
of the total angular momentum ^ASi + Si) evaluated in 
the local frame. Given that J > 0, in the ground state, 
Si and Si are then parallel. 

This approximation is manifestly valid in the half 
metal limit when J 3> t. The ground state is then a 
mixture of states in which all sites are either singly occu- 
pied by an electron or unoccupied. The singly occupied 
sites with the maximum angular momentum 5 + (1/2) 
have the lowest energy while other states, and those with 
two electrons per site, have an energy which is higher 
by ~ J and hence have negligible weight in the ground 
state. The Wigner-Eckart theorem then dictates that all 
the matrix elements of s*j are equal to those of (5^/25). 

However, since spatially the magnetization rotates 
slowly, a much weaker inequality suffices. The adia- 
batic theorem simply demands that the transverse field 
Atjj = it(V9i ■ Tij)a y be small compared to the longitu- 
dinal field Jm s . The wall rotates by 7r over a distance w 
so V9i ~ ir/w and At^ ~ iirt{a/w), and required is, 

nt(a/w) <C Jm s (5) 

which since, e.g., for Permalloy w/a ~ 10 3 is typically 
well satisfied. The conduction electron magnetization 
comprises two components with, by definition, the (mi- 
nority) majority conduction electrons (anti-) parallel to 
the axis of quantization, i.e., the direction of the local 
spin. In the local frame the majority (minority) elec- 
trons have cr 2 = +1 (a z = —1), so that it follows that 
when the adiabatic theorem is satisfied, 

Si = a z (Si/2S), (6) 

independent of the details of the electronic structure, etc. 

The non-interacting approximation is therefore to in- 
deed ignore At^. If it was not for the cosine in tij = 
tcos(V9i ■ fij) the problem is then identical to that with- 
out the wall. Using again (V0j • fij) ~ Tta/w, the correc- 
tion is ~ tir 2 (a/w) 2 which with (a/w) ~ 1CP 3 might be 
safely ignored, i.e., at this level the electronic structure 
is unchanged by the wall and since the reflection prob- 
ability is small, there is a negligible pressure exerted by 
the conduction electrons. However, in the spin sector, 



3 



this correction is important. Evaluating the coefficient 
of (V0j • Tij) 2 in first order perturbation theory leads to, 
the renormalization J s = J° + (x't/2S 2 ), of the exchange 
coupling (liij. The effective concentration x' — (c|c-). 

At this lowest level of approximation the wall is sta- 
tionary and is put in motion only when the torque trans- 
fer Ai^-term is accounted for as a perturbation. Even 
with this term, the effective Hamiltonian, Eqn. (J2J) , is ev- 
idently of single particle nature. Consider first the single 
particle description of the majority spin electrons. In or- 
der to account for the torque transfer term proportional 
to Ate, use is made of Eqn. |(BJ. To this end, it is noted, 

for majority electrons J^aa' C L°W C .7<7' = K s i c n c j-\ ~ 



c iT c J'T s j 



Then by virtue of Eqn. |JBJ, e.g., s~l 



(2S)- 1 /\, and E,^ cla yaa ,c ja , = i{2S)- l / 2 c\^{b\- 
6 ). Combining this with the similar result for the mi- 
nority electrons, the Ai^-term becomes, 



t 



2(2S) 1 / 2 



<T(V6en j )cl c jAbi-b i j )+H.c., (7) 



which couples the spin current to the magnons and re- 
flects the entire torque transfer process. The necessary 
correction is obtained, directly, by taking the expectation 
value with respect to the conduction electrons, whence 
Eqn. |(7J| reduces to 



. Hj s a 2 
' 2eS 



^(vdi-f^sy/w-bD, (8) 



where j s is the spin current. Here b i+1 is replaced by b\, 
an approximation valid in the continuum limit w 3> a. 
The quantity (2S) 1 ^ 2 (b i — b\) ~ Si V and Tt T corresponds 
to an effective field which is strictly transverse. The ap- 
pearance of such a term linear in (b i — b\) signals that 
there is no time independent solution. 

Comparing H T with the definition of p z makes evident 
that Ti. T oc j s p Zl i.e., the current only couples to the 
collective coordinate of the wall via its momentum. The 
effect of TL T is then to put the wall in motion along the z- 
direction, and it is necessary to study the time dependent 
Schrodinger equation: ih(d/dt)ip(ri,t) = Ti,ip(fi,t). The 
effect of adding the rotations, r, which displace the wall 
is ih{d/dt) -> ihrid/dty- 1 = ih(d/dt) - h{d9i/dt)M iy , 
where Mi = Si + Si is the total spin angular momentum. 
This generates a second purely transverse field term in 
the effective Hamiltonian: 



S 



using the fact that (Mi/M) = (Si/S), i.e., that all mag- 
netizations are parallel. Thus when 



dt 



v - 



dz 



J s a J 
2M 



(9) 



the effective fields generated by the spatial and temporal 
rotations of the axes of quantization cancel each other. 
Given 6> t = 0(z t ) = 2 cot" 1 e"^*/^ for j s = 0, the new 
ground state has 9i = 6{zi — Vot), i.e., the wall moves 
without distortion and without tilting or twisting. It 
is easy to show that the result vq = (j s v c /2M) reflects 
the conservation of the z-component of the total angular 
momentum, i.e., that the net spin current, 2j s carried 
towards the wall by the electrons equals the change in 
the angular momentum, js — Mv/v c , of the wall due to 
its motion. Since the conduction electrons are polarized 
j s = pj is related to the charge current j by some material 
determined parameter p. 

It is to be observed that by making these specific 
time dependent rotations the net transverse magnetic 
field has vanished from the problem. Thus the condi- 
tions for equilibrium in rotating frame are identical to 
those for the static problem, when j s — 0, and in par- 
ticular dV/d(f>i = since this is an equilibrium condi- 
tion. In order to verify this point of divergence with 
Bazaliy et al.0, it is useful to calculate this derivative 
directly. This requires the full SU(2) transformations. 
After some algebra the result is, Eqn. 10} with, A = 
(1/2) (V0iK'+ (1/2) sin 9i (V<kW+ (1 - cos 0i)(V&)o£ 
where the a' are defined in the local frame of reference. 
Directly, both new terms are zero for a simple domain 
wall for which = 0, and dV/dfc = o[y|. 

In terms of the Landau-Lifshitz equations, in the lab- 
oratory frame, 

^^^i-0 s . V)M = gp B M x B (10) 
where the effective field is (K± = for simplicity): 

B = ^r = JS 2 a 2 V 2 M+^{z-M)z. (11) 
3M M 2 

This is formally similar to the result obtained by Baza- 
liy et al.|9j, in the half metal limit, except for the fact 
that the term proportional to j s does not arise from a 
finite dV/d(j>, via Eqn. I|ll|) . Also this torque transfer 
term arises from At^j which is simply not included in 
that work. Equation (fTUf) defines the "paticle derivative" 
DM/Dt, i.e., that at a fixed position in the moving wall. 

This different origin of the torque transfer term is of 
particular importance when relaxation is accounted for. 
With j s = 0, relaxation is traditionally included through 
a Gilbert term —{a/M)M x (dM/dt), parameterized by 
a, and which assures the system relaxes until M is paral- 
lel to the internal field B, and which corresponds to an ab- 
solute minimum of V{9i,4>i}. This is an evident require- 
ment of the second law of thermodynamics. In order that 
this law be satisfied when j s is finite, the Gilbert term 
must involve rather DM/Dt and the Landau-Lifshitz- 
Gilbert equations are, 



DM 
~Dt 



gpsM x B 



a - DM 
— M x . 

M Dt 



(12) 
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with the present Eqn. Ijlljl. 

The justification, at the microscopic level, of this form 
of relaxation is complicated. For the micron sized sys- 
tem it is widely assumed that the principal loss mecha- 
nism corresponds spin-orbit scattering of the conduction 
electrons 0- Above a few tens of degrees, the Yafet-Elliot 
|l5|. spin-orbit scattering off phonons is almost certainly 
dominant. However, here there is little interest in ob- 
taining a detailed expression for (1/r) and so rather con- 
sidered is the simplest model, 

H so = *A]T££e^-^<w • (k x k')cy %ltT , 

' kk< 

which represents such scattering off impurities, located 
are r*i, and which enters the theory in the same man- 
ner as does phonon scattering. In fact, if the concen- 
tration c of "impurities" is taken to be ~ T/Td where 
Td is the Deby temperature, this mimics the effect of 
phonons at modest temperatures. The problem is treated 
in the rotating frame since this reduces to zero the gauge 
fields created by the current j s . Now in H so it is nec- 
essary to make the replacement a — > rar^ 1 . The an- 
gular speed of rotation, ~ irv /w ~ 10 3 rads/s, is by far 
the lowest frequency in the problem and hence it can be 
assumed that, u>o — @i is a constant during the relax- 
ation process. Further, in practice, the wall width w is 
large compared to the mean free path, and so the spa- 
tial dependence of ujq can be ignored. The problem is 
then very similar to the linear response to an applied ra- 
dio frequency field. In the context of dilute magnetic 
alloys, this has been studied in some detail[lfij. The 
Gilbert term is found to be of the correct form with, 
in the present notation, a = \/{uj s t) where uj s — g^BB 
and, (1/r) = Wf_>i — Wj, _►•[•, and where, e.g., Wf_>i is 
the transition rate for conduction electrons to flip from 
up to down spin given essentially by the "Golden Rule" . 
Assuming a random concentration of "impurities", e.g., 
= c{-K\*/h)Y, m ,[{kzK - k x k' z )H{e % - e p ) + 
(l/A)((k y k' z - k z k' y f + (k x k' y - k y k' x f)8{e % - eg, -Huo)) + 
(l/4)((k y k' z - k z k' y f + (k x k' y - k y k' x f)5{^ k - e- k , + fiL*))], 
where eg is the energy of an electron with wave vector k. 
The only effect occasioned by the passage to the rotat- 
ing frame is the appearance of hu>o in the delta functions 
with quantum corrections to the form of the Gilbert term 
which are negligible since hu>o <C ksT. The modifications 
to the expression for (1/r) are also negligible since, even 
for the true phononic mechanism, Hujq is much smaller 
than the relevant energy scale A^T. 
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